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(57) Abstract 

A method for simulating fluid flow within a three-dimensional object having first and second generally opposed surfaces including 
matching each element (DF) of said first surface with an element (CE) of said second surface between which a reasonable thickness may 
be defined, wherein matched elements of said first surface constitue a first set of matched elements and matched elements of said second 
surface constitute a second set of matched elements specifying a fluid injection point performing a flow analysis using each set of said 
matched elements, whereby said injection point is linked to all locations on said first and second surfaces from which flow may emanate 
such that resulting flow fronts along said first and second surfaces are synchronized. 
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METHOD FOR MODELLING THREE DIMENSIONAL OBJECTS 

AND SIMULATION OF FLUID FLOW 

This invention relates to a method for modelling solid 
5 objects, particularly for use in the simulation of fluid 
flow, to be used for example to simulate prototypes before 
production. In a preferred embodiment the method is used 
in the design of articles to be manufactured by injection 
molding, preferably from molten plastic materials. 

10 

The modelling of solid objects is employed in various 
fields. Such modelling is used, for example, in the 
simulation of injection molding. In that field, it is 
widely recognized that the filling and packing phases of \, 

15 injection molding have a significant effect on the visual 
and mechanical properties of a molded object. Simulation 
is employed to analyse proposed shape b and injection 
points, and thus the final quality of the ultimate article. 
A requirement of any injection mold 1b that it can be 

20 filled with molten polymer given the pressure limits of a 
real injection molding machine . Simulation can provide 
information as to whether the mold can be filled and the 
fill pattern that will be achieved. By using simulation, it 
is possible to determine optimum gate locations and 

25 processing conditions. It is possible to predict the 

location of weld lines and air traps . Economic benefit is 
derived from simulation because problems can be predicted 
and solutions tested prior to the actual creation of the 
mold. This eliminates costly re-working and decreases the 

30 time required to get an object into production. 

Simulation technology has been developed and generally uses 
finite element /finite difference techniques to solve the 
governing equations of fluid flow and heat transfer. In 
35 order to minimize the time required for analysis and hence 
the required computer resources, the He le- Shaw 
approximation is invoked to simplify the governing 
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equations . It has bees found that tills simplification 
provides sufficient accuracy for injection molding but does 
create the need for specific modelling of the computational 
domain. 

5 

Injection molding is an excellent process for repetitively 
manufacturing large numbers of objects or parts having 
complicated geometries. A characteristic of injection 
molded components id that the thickness of the wall is 
10 generally a small fraction of the overall length, of the 
component* In view of the low thermal conductivity of 
plastics, this physical characteristic is essential to 
achieve the rapid cycle times that make the process so 
attractive. 

15 

The flow of melt in an injection mold is determined by the 
familiar conservation laws of fluid mechanics. Solution of 
the equations in their full generality presents several 
practical problems. Owing to the characteristically thin 

20 walls of molded components, however, it is possible to make 
some reasonable assumptions that lead to a simplif ication 
of the governing equations. These simplified equations 
describe what is called Hele-Shaw flow and may he readily 
solved in complex geometries using an appropriate numerical 

25 technique such as the finite element and/or finite 
difference method. 

Injection molding simulation is now routinely regarded as a 
desirable aspect of plastic part design. Similarly, 

30 improved computer aided drafting (CAP) technology has led 
to the widespread use of surface and solid modelling. 
Associated advantages of this are the ability to better 
visualize an object, to use numerical cutting, and the 
possibility of achieving more concurrency in engineering 

35 design and manufacture. When using the Hele-Shaw 

approximation, plastic CAB analysis still requires the use 
of a surface model, representing the midplane of the real 
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component, which is then meshed with triangular or 
quadrilateral elements to which suitable thicknesses are 
ascribed. The preparation of such a mesh can take a 
considerable amount of time, and requires substantial user 
5 input; owing to the labour intensive nature of this step, 
model preparation requires the greatest share of time when 
performing a molding simulation and makes this technique 
time consuming. In addition, as model preparation is an 
interactive task, it has a higher cost associated with it 
10 than simply running a computer program. 

This conventional technique for the simulation of three 

dimensional bodies is illustrated in figures 1 to 3 of the 

accompanying drawings in which: 
15 Figure 1 is a representative example of an 

article to be manufactured by injection molding in the form 

of a T-shaped component; 

Figure 2 is a midplane representation of the 

component of figure 1 that is meshed; and 
20 Figure 3 is a view of the major results from a 

conventional analysis showing the position of the advancing 

melt at given times. 

One solution to the above shortcomings is to avoid the use 
25 of the He le- Shaw equations and solve the governing 

equations in their full generality. This has inherent 
problems owing to the thin walled nature of injection 
molded objects and parts. To perform such an analysis, the 
region representing the mold cavity into which molten 
30 polymer will be injected must be divided into small sub- 
domains called elements. Usually these elements are of 
tetrahedral or hexahedral shape. This process of sub- 
division is called meshing and the resultant network of 
tetrahedra or hexahedra the mesh. Owing to the complicated 
35 shape of many injection molded objects and parts it is 
generally not possible to automatically mesh the cavity 
with hexahedral elements. It is possible, however, to mesh 
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the domain automatically with tetrahedral elements. The 
thin walled nature of injection molded objects and parts 
means that the plastic is subject to a huge thermal 
gradient in the thickness direction of the component. This 
5 requires that there be a reasonable number of elements 

through the thickness. Using existing meshing technology, 
the result is a mesh consisting of hundreds of thousands or 
even millions of elements. The high number of elements 
makes the problem intractable for any but the fastest super 
10 computers. These are rarely found in industry, being 

extremely costly to purchase and maintain. Thus, although 
three dimensional simulation provides a solution that 
avoids the requirement of a midplane model, it is not as 
yet a practical solution. 

15 

it is an object of the present invention, therefore, to 
provide a method for the simulation of flow in a three 
dimensional object that can produce simulations 
substantially automatically, without requiring the solution 
20 of the governing equations in their full generality. 

According to a first broad aspect of the present invention, 
therefore, there is provided a method for simulating fluid 
flow within a three dimensional object having first and 

25 second generally opposed surfaces including: 

matching each element of said first surface with 
an element of said second surface between which a 
reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 

30 matched elements and matched elements of said second 
surface constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using each set of said 
matched elements, whereby said injection point is linked to 

35 all locations on said first and second surfaces from which 
flow may emanate such that resulting flow fronts along said 
first and second surfaces are synchronized. 
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It should be noted that, in this context, generally opposed 
surfaces may be parallel, inclined to each other, meet at 
an acute angle or otherwise, and need not be planar, 
5 provided that a reasonable thickness or thicknesses may be 
assigned to the space between the surfaces. 

Preferably said injection point is first linked to all said 
locations from substantially the commenc emen t of said flow 
10 analysis* 

Preferably said injection point remains so linked at 
substantially all times in said flow analysis subsequent to 
said first linking. 

15 

Thus, although the injection point must be linked to all 
these locations at some point in the flow analysis, and 
essentially continuously thereafter, the injection point 
nay remain unlinked initially. Although this may decrease 
20 the accuracy of the final analysis adversely, it may 

nevertheless allow the analysis to yield useful results. 

Preferably said injection is one of a plurality of 
injection points. 

25 

Preferably said synchronization of said flow fronts is 
checked periodically. 

Preferably said flow fronts are synchronized if found not 
30 to be or no longer synchronized. 

Preferably said checking is performed at each time step. 

Thus, although for simple objects it is sufficient to 
35 synchronize flow fronts by linking the injection point or 
points to all locations, more complex parts require this 
synchronization at each calculation&l time step. 
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Preferably said first and second generally opposed surfaces 
are one of a plurality of pairs of generally opposed 
surfaces. 

5 

Thus, where the existing techniques utilize a single, 
mi dp lane representation of the object in which flow is to 
be modelled, and perform a simulation with this 
representation, the method according to the present 

10 invention utilizes only the outer surfaces defining the 

three dimensional object to create a computational domain. 
These correspond to representations (in a preferred 
embodiment, meshed representations) of the domain in which 
flow is to be simulated, and would comprise, for example, 

15 meshed representations of the top and bottom surfaces of a 
part. Thus, in such an example the invention could be said 
to utilize an outer skin mesh rather than a midplane mesh. 
Elements of the two surfaces are matched, based on the 
ability to identify a thickness between such elements. An 

20 analysis, substantially along conventional lines (by means, 
for example, of the Hele-Shaw equations) , is then performed 
of the flow in each of these domains in which flow is to be 
simulated, but linked to ensure fidelity with the physical 
reality being modelled. 

25 

Preferably any unmatched elements of said first and second 
surfaces, being elements that could not be matched, are 
assigned thicknesses being the average of the thicknesses 
of adjacent matched elements where such adjacent matched 
30 elements exist, or of adjacent unmatched elements where 
such adjacent matched elements do not exist and said 
adjacent unmatched elements have been assigned thicknesses. 

Thus, eventually all matched and unmatched elements of the 
35 first and second surfaces may be assigned thicknesses. 

Preferably each element of an edge surface, being a surface 
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between said first and second surfaces, and adjacent to 
either of said first or second surface is assigned a 
thickness proportional to the thickness of the element of 
said first or second surface to which said element of said 
5 edge surface is adjacent. 

Preferably each said element of an edge surface is assigned 
a thickness between 0.5 and 1.5 times said thickness of the 
element of said first or second surface to which said 
10 element of said edge surface is adjacent. 

Preferably each said element of an edge surface is assigned 
a thickness between 0.7 and 0.9 times said thickness of the 
element of said first or second surface to which said 
15 element of said edge surface is adjacent . 

Still more preferably each said element of an edge surface 
is assigned a thickness 0.75 times said thickness of the 
element of said first or second surface to which said 
20 element of said edge surface is adjacent. 

Preferably each element of an edge surface not adjacent to 
said first or second surface is assigned a thickness being 
the average of the thicknesses of adjacent elements of said 
25 edge Burface that have been assigned thicknesses. 

Assigning a thickness/ therefore, to every element of both 
surfaces and of edge surfaces ensures that the greatest 
accuracy will be achieved in the simulation. 

30 

Preferably flow is simulated at a rate directly 
proportional to a desired flow rate for the object. 



35 



The desired flow rate will usually be a flow rate selected 
by a user based simply on the volume of the object to be 
filled and the desired fill-time. 
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Preferably eaid rate is proportional to the ratio of 
computational domain volume of said object to real volume 
of said object. 

5 More preferably said rate is substantially equal to the 
ratio of said computational domain volume to said real 
volume. 

This approach may be used to compensate for the 
10 modification of the effective volume of the object being 
modelled/ owing to the use of two domains (for example, 
associated with the top and bottom surfaces), rather than 
one meshed domain as in the midplane approach. The 
computational domain volume is thus this effective volume 
15 used for the purposes of modelling the object, which most 
preferably will be approximately double the real volume. 

Preferably the method is performed with first and second 
representations of said first and second surfaces 
20 respectively comprising first and second meshes or lattices 
respectively, wherein said elements are interstices of said 
first and second meshes or lattices. 

Preferably said elements are triangular or quadrilateral. 

25 

Preferably, when said elements are triangular, said 
elements are substantially equilateral. 

Preferably said method includes creating said first and 
30 second representations. 

In one embodiment the method includes creating improved 
representations of said first and second surfaces, whereby 
said elements are elements of said improved representations 
35 and said method is performed with said improved 
representations . 
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Preferably said first and second representations are, or 
are a part of, a representation or representations for 
stereolithography of said object • 

5 Thus, computer representations of the two surfaces may be 
provided as input to the method, may be created by the 
method, or - if provided with representations of the 
surfaces - the method may create improved representations. 
As described above, particularly preferred representations 
10 include those with elements comprising equilateral 
triangles . 

in one embodiment of the invention the method described 
above is performed by a computer running a computer program 
15 encoding said method for simulating fluid flow. 

Preferably said method optionally includes corrections for 
non- isothermal temperature fields and/ or non-Newtonian 
fluids. 

20 

Thus, known techniques for including the effects of non- 
isothermal temperature fields and non-Newtonian fluids may 
be included so that more physically faithful simulations 
may be performed, as well as faster simulations without 
25 these corrections when speed is desirable, even at the 
expense of precision* 

According to a second broad aspect of the present invention 

there is provided a method for simulating fluid flow within 
30 a three dimensional object having first and second 

generally opposed surfaces including: 

providing or creating first and second 

representations of said first and second surfaces 

respectively, 
35 creating first and second improved 

representations from said first and second representations 

respectively. 
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matching each element of paid first improved 
representation of said first surface with an element of 
said second improved representation of said second surface 
between which a reasonable thickness may be defined, 
5 wherein matched elements of said first improved 

representation constitute a first set of matched elements 
and matched elements of said second improved representation 
constitute a second set of matched elements, 
specifying a fluid injection point, 

10 performing a flow analysis using each set of said 

matched elements, whereby said injection point is linked to 
all locations on said first and second improved 
representations from which flow may emanate such that 
resulting flow fronts along said first and second improved 

15 representations are synchronized. 

Preferably the first and second representations are, or are 
a part of, a representation or representations for 
stereo lithography of said object. 

20 

Preferably the first and second improved representations 
comprise small equilateral triangular elements. 

According to a third broad aspect of the invention there is 
25 provided a computer provided with or running a computer 
program encoding the method for simulating fluid flow 
described above. 

In a fourth broad aspect of the invention there is provided 
30 a computer storage medium provided with a computer program 
embodying the method for simulating fluid flow described 
above . 

According to a fifth broad aspect of the present invention 
35 there is provided a method for modelling a three 
dimensional object including: 

specifying first and second generally opposed 
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surfaces of said object, 

forming first and second representations of said 
first and second surfaces respectively, wherein said first 
and second representations each comprise a plurality of 
5 elements, 

matching pairs of elements of said first and 
second surfaces between which a reasonable thickness may be 
defined. 

10 Preferably the first and second representations comprise 
first and second meshes or lattices respectively, wherein 
said elements are interstices of said first and second 
meshes or lattices. 

15 Thus, the elements of the first representation are 

interstices of the first mesh, and the elements of the 
second representation are interstices of the second mesh. 

Preferably said elements are triangular or quadrilateral. 

20 

Preferably, when said elements are triangular, said 
ele m ents are substantially equilateral. 

Preferably each element of each of said matched pairs of 
25 elements is assigned respectively said thickness. 

Preferably any unmatched elements of said first and second 
surfaces, being elements that could not be matched, are 
assigned thicknesses being the average of the thicknesses 
30 of adjacent matched elements where such adjacent matched 
elements exist, or of adjacent unmatched elements where 
such adjacent matched elements do not exist and said 
adjacent unmatched elements have been assigned thicknesses. 

35 Preferably each element of an edge surface, being a surface 
between said first and second surfaces, and adjacent to 
either of said first or second surface is assigned a 
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thickness proportional to the thickness of the element of 
said first or second surface to which said element of said 
edge surface is adjacent. 

5 Preferably each said element of an edge surface is assigned 
a thickness between 0.5 and 1.5 times said thickness of the 
element of said first or second surface to which said 
element of said edge surface is adjacent. 

10 Preferably each said element of an edge surface is assigned 
a thickness between 0.7 and 0.9 times said thickness of the 
element of said first or second surface to which said 
element of said edge surface is adjacent. 

15 Still more preferably each said element of an edge surface 
is assigned a thickness 0.75 times said thickness of the 
element of said first or second surface to which said 
element of said edge surface is adjacent. 

20 Preferably each element of an edge surface not adjacent to 
said first or second surface is assigned a thickness being 
the average of the thicknesses of adjacent elements of said 
edge surface that have been assigned thicknesses. 

25 According to a sixth broad aspect of the present invention 
there is provided a method for simulating fluid flow within 
a three dimensional object having first and second 
generally opposed surfaces including: 

matching each element of said first surface with 

30 an element of said second surface between which a 

reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of Baid second 
surface constitute a second set of matched elements , 

35 specifying a fluid injection point, 

performing a flow analysis using each set of said 
matched elements, and 
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synchronizing flow fronts resulting from said 
flow analysis along said first and second surfaces. 

Preferably said flow fronts are synchronized from 
substantially the commencement of said flow analysis. 

Alternatively said flow fronts are first synchronized after 
the commencement of said flow analysis. 

According to. a further broad aspect of the present 
invention there is provided a method for simulating fluid 
flow within a three dimensional object having first and 
second generally opposed surfaces including: 

matching each element of said first surface with, 
an element of said second surface between which a 
reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements, 

specifying a fluid injection point, 

performing a flow analysis using said first set 
of matched elements, 

adapting said flow analysis to said second set of 
matched elements, and 
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synchronizing flew fronts resulting from said 
flow analysis and said adaptation of said flow analysis 
along said first and second surfaces . 

5 Thus, it is possible in some cases to perform the method 
using only one of the sets of elements, and mapping the 
result onto the other set, and performing minor adaptation 
of the analysis to allow for any differences between the 
two surfaces. 

10 

Preferably the method is performed with first and second 
representations of said first and second surfaces 
respectively comprising first and second meshes or lattices 
respectively, wherein said elements are interstices of said 
15 first and second meshes or lattices. 

Preferably the elements are triangular, and most preferably 
equilateral triangles. 

20 Preferably the elements are quadrilateral. 

Preferably the method includes creating said first and 
second representations, and more preferably creating 
improved representations, either from the original surfaces 
25 or from the first and second representations, and 

performing the method with elements of the representations 
or improved representations. 

In all the above methods according to the present invention 
30 aspects it is preferred that the synchronization comprises 
matching pressure and temperature. 

According to particular aspect of the present invention 
there is provided a method for simulating fluid flow within 
35 a three dimensional object having first and second 
generally opposed surfaces including: 

matching each element of said first surface with 
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an element of said second surface between which a 
reasonable thickness may be defined, wherein matched 
elements of said first surface Constitute a first set of 
matched elements and matched elements of said second 
5 surface constitute a second set of matched elements, 

specifying a plurality of fluid injection points, 
performing a flow analysis using each set of said 
matched elements, whereby said injection points are linked 
to all locations on said first and second surfaces from 
10 which flow may emanate such that resulting flow fronts 
along said first and second surfaces have matched flow 
rates. 

A preferred embodiment of the invention will be described, 
15 by way of example, with reference to the accompanying 
drawings, in which: 

Figure 4 is an example of a stereolithograpby 
mesh of the T-shaped object of figure 1; 

Figure 5a shows an initial mesh used in a 
20 remeshing method according to the present invention; 

Figures 5b- f dhow the mesh of figure 5a 
progressively remeshed according to the remeshing method; 

Figure 6 shows the mesh of figure 4, remeshed 
according to the preferred embodiment of the present 
25 invention; 

Figure 7 a is a view of a mesh produced by 
conventional modelling techniques for a flat plate; 

Figure 7b is a view of a corresponding mesh 
produced for a flat plate according to a preferred 
30 embodiment of the present invention; 

Figure 8 illustrates surface matching for a flat 
plate according to the preferred embodiment; 

Figure 9 is a view similar to figure 8 
illustrating in cross-section the edges, matched surfaces 
35 and unmatched surfaces of another T-shaped object according 
to the preferred embodiment; 

Figure 10 is a cross -sectional view illustrating 
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the edges, matched surfaces, matched tapered surfaces, 
unmatched surfaces and unmatched edges of yet another T- 
shaped object according to the preferred embodiment; 

Figure 11 is a view of a pimple object showing 
5 injection points; 

Figure 12 is a view of an object comprising 
intersecting plates with a user selected injection point 
according to the preferred embodiment; 

Figure 13 is a enlarged view of the object of 
10 figure 12 showing the injection points required to allow 

flow to emanate correctly from the selected point according 
to the preferred embodiment; 

Figure 14 is a cross-section of a plate, with a 
flow front advancing from the left and leading a little on 
15 the top surface according to the preferred embodiment; 

Figure 15 illustrates in a cross-sectional view 
the synchronization of flow in a ribbed plate, with the 
flow splitting and also filling a vertical rib according to 
the preferred embodiment; 
20 Figure 16 illustrates in a cross -sectional view 

the assignment of identical pressures at nodes to 
synchronize flow fronts at a rib according to the preferred 
embodiment ; and 

Figure 17 is a view similar to figure 3 generated 
25 according to the preferred embodiment of the present 
invention. 

A method for simulating fluid flow within a three 
dimensional object according to the present invention 
30 involves three main steps: 

1) Preparation of geometry, 

2) Selection of injection points, and 

3) Analysis. 

Before these steps are effected, however, the input for the 
35 simulation is prepared. Mesh generators for producing 

finite element meshes are generally an expensive addition 
to the base CAB system. Consequently many CAB 
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installations have no finite element analysis (PBA) mesh 
generation capability* xn order to facilitate the 
implementation of the invention, the method has been 
designed to avoid the need for a xnesb generator. This is 
5 done by using the type of mesh available for 

stereo lithography . Stereolithography, a prototyping 
technology that is widely employed in industry, is used to 
produce prototype components for assembly or evaluation. 
It requires a stereolithography apparatus and takes as 

10 input a computer file known as a "SLA" file. The ability 
to output SLA files is common in CAP systems. However, the 
type of triangulation created is not good enough for FEA. 
Characteristic of the SLA file is the inclusion of several 
very long, narrow triangles of high aspect ratio as shown 

15 in figure 4* 

A method in accordance with the present invention includes, 
therefore, a remeshing algorithm (or remesher) that uses as 
input the triangles from a stereolithography SLA (or 
20 similar) file and remeshes to produce a mesh that may be 
Used in analysis. 

There are several steps involved in the input of a 
stereolithography mesh and its refinement into elements for 
25 analysis. These steps are: 

• Input of mesh, 

• Checking of mesh, 

• Classification of mesh into surfaces, and 

• Remeshing. 

30 Each of these is described below. 

Input of Mesh 

In the preferred embodiment, the input mesh is of the form 
used for stereolithography. Alternatively, however, the 
35 input could be IOES surfaces (which would then be meshed 
internally), a surface mesh of quadrilaterals or a mix of 
quadrilaterals and triangles, tetrahedral or hexahedral 
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meshes • 

Checking of mesh 

The input mesh is checked to ensure that it is closed and 
5 orientable* The mesh is consistently oriented so that a 
unit normal to each element points outward from the 
internal volume defining the solid region that represents 
the object or part being modelled. The volume of this 
region and the surface area of the mesh are calculated and 
10 stored. The number of parts defined by the input mesh is 
also determined. 

V 

Classification of mesh into surfaces 

The following properties for each element in the mesh are 
15 noted or calculated, and stored: 

• Area/ 

• Normal at the element cent ro id, 

• Edge lengths of the element, 

• Internal angles at the vertices of the elements, 
20 - Adjacent elements (the elements attached to each 

side of the element under consideration) , and 

• Bending angle between adjacent elements (the 
angle between normals of adjacent elements). 

Nodal properties are then noted, calculated and stored. 
25 These are: 

• Measure of minimum curvature (calculated by 
looking at the elements attached to a node and noting the 
angle between each adjacent pair of elements, whereby the 
measure of minimum curvature is then the smallest such 

30 angle between attached pairs of elements), 

• Connecting elements (the set of elements 
connected to the node), and 

• Number of edges connected to the node. 

Surface edges are then calculated using a value of bending 
35 angle to group elements into surfaces. Surfaces so formed 
are then classified into 

• planar surfaces, and 
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• low curvature surfaces* 

Further classification is then done to produce high 
curvature surfaces and edges. The following properties of 
edges are noted or calculated and stored: 
5 • Length, 

• Bending angle , 

• Direction of bending (in or out), and 

• Adjacent Elements. 

Surface loops are then created. These are defined to be 
10 the oriented edges of surfaces. It should be noted that a 
surf ace with a hole in it will have two associated loops - 
one for the outer edges and one describing the interior 
hole. Loop properties are then noted or calculated and 
stored: 
15 • Length, and 

• Edges connected to loop. 

The following surface properties are noted or calculated 
and stored: 

• Perimeter, 
20 • Area, 

• Nodes in surface, 

• Elements in surface, 

• Edges, and 

» Minimum measure of curvature. 

25 

Remeshing 

With the surfaces classified, the mesh associated with each 
surface is refined using a bisection algorithm which 
bisects the longest side of a triangle and creates extra 

30 elements. This is illustrated in figures 5a~f . Figure 5a 
shows the initial mesh. The mesh is refined by defining a 
node at the mid point of the longest; element side and 
extending lines to vertices to define extra triangles. In 
figures 5b to 5f , the midpoint of the longest side (denoted 

35 by o) and the dotted lines extending from this point define 
the new elements. Remeshing continues until the elements 
satisfy a criterion on size. 
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The result of the remeshing algorithm when applied to the 
mesh shown in figure 4 is shown in figure 6. In this 
regard the ideal triangle is equilateral: it may be seen 
5 that the triangles in figure 6 are far more uniform and 
closer to this ideal. 

The method will thus accept as input stereolithography 
files (usually generated with a " .stl" filename extension) 
10 and other files that are similar such as "render" files 
(produced by "Pro -Engineer" with a " .sip" extension). 

In practice any file format that describes a mesh (of, for 
example, a lattice or quadrilaterals) covering the outer 
15 surfaces of the solid region may be used as input to the 

invention. The remesher is always used to improve the mesh 
quality and produce the triangular mesh for subsequent 
processing. 

20 The first step is then the preparation of the geometry, in 
which the geometry of the three dimensional CAD solid model 
is processed by conventional techniques to provide a 
surface mesh covering the outer surfaces of a solid body, 
rather than a mi dp lane mesh as required by the conventional 

25 simulation approach. This distinction is depicted in 

figures 7a and 7b. Figure 7a shows the midplane mesh for a 
simple plate required by the conventional technique, 
whereas figure 7b shows the mesh employed according to the 
present invention for the same object. Further, the flow 

30 front must be synchronized on both sides of the surface in 
which the material flows. If material is injected at the 
centre of the plate, then, in the conventional case, the 
material flows to the boundary of the part or object and 
stops without any difficulty. For the surface mesh on the 

35 solid model (figure 7b) the material will flow to the outer 
edge, across the edge and then under the top surface. This 
is clearly not indicative of the physical phenomenon of 
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flow in a plate and 00, according to the invention, there 
are imposed some constraints on the flow when using the 
surface mesh. These constraints are that the injection 
point must link to all surfaces from which flow will 
5 emanate and the flow front must be synchronized along the 
top and bottom meshes of the surface mesb. Preferably this 
link is established from the commencement of the flow 
analysis, but a successful (if soaewbat less accurate) 
analysis of the flow can still be performed is this linking 

10 commences subsequently. To perform these tasks the surface 
mesh must be categorized into surfaces wbich are then 
further classified to facilitate injection node selection 
and synchronization. Thus, according to the present 
invention, the surface mesh is then grouped into individual 

15 surfaces and classified. The categories of surfaces are: 

• Matched, 

• Unmatched, and 

• Edges • 

Hatched surfaces are defined to be tbose that have a 
20 relationship with another surface such that a notion of 
thickness between them can be sensibly defined. Figure 8 
illustrates this idea. Surfaces abed and efgh may clearly 
be associated so that the thickness t between them is 
defined. Surfaces abfe, begf , cdhg and daeh are edge 
25 Surfaces and unmatched. 

Figure 9 further clarifies these concepts. This figure 
shows a cross-section through a filleted rib. Here ab, cd 
and gh are on edge surfaces. Lines aj and ed are matched 

30 to be, fg and hi are matched and the curved sections ij and 
ef are unmatched. It should be noted tbat it is not 
possible to sensibly define a thickness of unmatched 
surfaces in the manner described above. As a final 
example, figure 10 illustrates a tapered rib cross -section. 

35 Here ab, cd and gh are edges as in figure 9. Lines aj and 
ed are matched to line be, and the curved lines ef and ij 
are unmatched as in figure 9. However fg and hi are still 
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considered matched despite the taper. Clearly, if the 
taper is extreme it is possible that the surfaces forming 
the tapered ribs may not be matched. 

5 When surfaces have been classified it is possible to assign 
to elements on matched surfaces a thickness equal to the 
distance between the surfaces. Elements on edge surfaces 
are given the thickness of the matched surfaces to which 
they are attached. Finally elements on unmatched surfaces 
10 are given the average thickness of surrounding elements on 
matched surfaces. 

According to the present invention, the mesh on each side 
of the object is analyzed. A key aspect in this is 

15 obtaining a similar flow front on each side of the matched 
surfaces. As the entry point for flow is the injection 
node and is selected by the user, the method ensures that, 
regardless of where the user chooses to inject plastic, all 
matched surfaces near the point are tied to the injection 

20 node. "Tied" implies that the same pressure and 

temperature of melt at the injection node is given to all 
other tied nodes. 

Figure 11 shows a cross-section of a rectangular plate. 

25 Point A is where the user chooses to inject plastic, point 
AA is the corresponding point on the other side of the 
surface. Point AA is automatically determined in 
accordance with the invention using the information from 
the classification of surfaces described above. Injection 

30 at this point ensures that the flow is matched on top and 
bottom surfaces. Similarly, point B is the point chosen by 
the user and SB is the corresponding point determined by 
the method in accordance with the invention to ensure that 
flow is matched along top and bottom surfaces. 

35 

The selection of injection node(s) can be quite complex. 
This is illustrated in figure 12, which shows an object 
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composed of intersecting plates. The small arrow near the 
point of intersection of the plates indicates the user 
selected injection point. Figure 13 is an enlarged view of 
the injection area and shows the actual injection points 
5 required to allow the flow to emanate correctly from the 
selected point. In this case eight injection points are 
required one of which is supplied by the user and the 
other seven are determined by the invention. 

10 This method identifies which points need to be linked to 
the selected injection point with the concept of matching. 
For a given injection point apecified by the user, a set 
consisting of all surfaces which are matched to all 
surfaces connected to the injection point is defined. The 

15 following steps are then performed: 

1* For each surface in this set, the point closest 

to the selected injection point is noted. These closest 
points form a set of points that are potentially to be 
linked to the user selected injection point. 

20 2. Bach member of the set of potential injection 

points is checked to see if it is already an injection 
point specified by the user. If it has already been so 
determined, it is omitted from the set. 

3. For each potential injection point remaining, a 

25 set consisting of all surfaces which are matched to all 
surfaces connected to the potential injection point is 
defined. 

These steps are repeated until the set of potential points 
30 is unchanged. The set then forms the set of points to be 
linked to the user specified injection point. 

This procedure may be readily understood for the simple 
geometry of two crossed plates, for which the user 
35 specifies injection at the intersection of the plates, in 
the event that the surfaces connected to the specified 
injection point (s) are not matched to any other surfaces. 
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no attempt is made to link other points to the injection 
point (e). Instead, flow is allowed to emanate until it 
reaches points attached to surfaces that are matched to 
opposing surfaces. The flow fronts are then synchronised 
5 on the opposing surfaces by assigning links between points 
on the opposing surfaces in a manner similar to that 
described above for linking points to the injection point. 

Analysis uses a He le- Shaw type solver. The injection flow 
10 rate is approximately double that for the volume of the 

object as each element in the surface mesh has a thickness 
equal to the gap thickness and so the mesh volume is 
approximately double that of the object. The time step for 
flow front advancement is determined by the geometric 
15 connectivity of the control volumes attached to each node. 

The solver has the ability to synchronize flow on matched 
surfaces during analysis. Figure 14 shows a cross-section 
of a plate where it is assumed that the flow front is 

20 advancing from the left and is leading a little on the top 
surface. By definition each node may be empty, partially 
filled or filled. On each time step, whether all the nodes 
belonging to an element are filled is checked. The instant 
this occurs the element is defined to be the current last 

25 element to fill and a line is drawn from its centroid to 
its matched element (the matched element is known from the 
surface classification described above) defining the check 
point. This "check line" is indicated in figure 14. The 
nearest node to the check point is then assigned the same 

30 pressure as the last node to fill, as indicated in figure 
14. Flow synchronization is done for all matched surfaces 
and all elements near the flow front (on top and bottom 
surfaces) are checked at each time step. 

35 Synchronisation is not performed on unmatched surfaces. 



An example highlighting the role of synchronization is 
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shown in figure 15, which shows the cross-section of a 
ribbed plate. The flow will split physically and fill the 
vertical rib as well as continuing to the right. Owing to 
the use of the surface mesh it is necessary to impose 
5 synchronization or the flow will continue up through node F 
to the top of the rib and then down to B and C. This is 
performed by using matching. Since the element comprising 
nodes D and F and the element comprising nodes C and E will 
be matched (as they are associated with matched surfaces) 

10 it is possible to assign the same pressures at nodes to 

synchronize the fronts. This is illustrated in Figure 16. 
Here a line from the centroid of the element comprising 
nodes D and F is constructed to its matched element 
(comprising nodes C and E) to define the check point. As 

15 the nearest node to the check point is C, it is given the 
same pressure as node r> for all subsequent time. Flow now 
emanates from node C as shown and so the flow goes up the 
rib on both sides as required. When establishing links for 
synchronization, links are only made to empty nodes to 

20 which are attached only empty nodes. 

A sample result of the analysis according to the present 
invention is shown in figure 17. This figure corresponds 
to figure 3, which shows the results of the comparable, 

25 conventional analysis of the same solid model shown in 
figure 1. As the present invention uses a surface mesh 
rather than a midplfcne mesh, the results are displayed on 
the solid model (figure 17) rather than on the midplane 
representation as in conventional analysis (figure 3). 

30 Apart from being more natural for the operator, the display 
is visually more meaningful. 

Fundamental quantities calculated by the invention, 
therefore , include : 
35 • the position of the melt front at any time (known 

as "filling contours"), 

* the pressure at the injection node and the 
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pressure distribution throughout the plastic at any time 
during filling or packing phases (known as a "pressure 
distribution"), and 

• the temperature of the plastic melt at any point 
5 within the part at any time during filling or packing 
phases (known as a "temperature distribution"). 
These may be displayed directly or processed to provide 
derived information to the user of the program, such as 
distributions of shear rate, shear stress, velocity, 

10 viscosity, direction of flow and orientation of reinforcing 
material. On the basis of this derived information and the 
fundamental calculated quantity data, changes to the 
component geometry, position of injection locations 
(gates), processing conditions used by the injection 

15 molding machine or material for molding may be evaluated 
for their efficacy in improving the equality or 
manufacturability of the part. This is generally an 
iterative process in which an analysis is performed, 
results are considered and changes made to bring about an 

20 improvement. Another analysis is then performed to ensure 
that the changes have indeed improved the results. 
Frequently changes to the part geometry are revealed. 
These are made in the CAD system and the revised model is 
subjected to further analysis. It is this aspect which is 

25 assisted by the invention as the need for the designer to 
recreate a new model for analysis every time a change is 
made is removed. 

Thus, an increase in the speed of evaluating of proposed 
30 part shapes is achieved. 

Modifications within the spirit and scope of the invention 
may readily be effected by persons skilled in the art. For 
example, linking may be performed in alternative methods 
35 such as, where multiple injection points are selected by a 
user, assigning the same flow rate to the linked nodes. It 
is to be understood, therefore, that this invention is not 
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THE CLAIMS DEFINING THE INVENTION ARB AS FOLLOWS: 

1. A method for simulating fluid flow within a three 
dimensional object having first and second generally 

5 opposed surfaces including: 

matching each element of said first surface with 
an element of said second surface between which a 
reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 
10 matched elements and matched elements of said second 
surface constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using each set of said 
matched elements, whereby said injection point is linked to 
15 all locations on said first and second surfaces from which 
flow may emanate such that resulting flow fronts along said 
first and second surf aces are synchronized ♦ 

2. A method as claimed in claim 1 wherein paid injection 
20 point is first linked to all said locations from 

substantially the commencement Of said flow analysis. 

3. A method as claimed in either claim 1 or 2 wherein said 
injection point remains linked to all said locations at 

25 substantially all times in said flow analysis subsequent to 
being first so linked* 

4 . A method as claimed in any one of the preceding claims 
wherein said injection point is one of a plurality of 

30 injection points. 

5 . A method as claimed in any one of the preceding claims 
wherein said synchronization of said flow fronts is checked 
periodically. 

35 

6. A method as claimed in claim 5 wherein said checking is 
performed at each time step. 
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7. A method as claimed in any one of the preceding claims 
wherein said flow fronts are synchronised if found not to 
be or no longer to be synchronized. 

5 

8. A method as claimed in any one of the preceding claims 
wherein said first and second generally opposed surfaces 
are one of a plurality of pairs of generally opposed 
surfaces . 

10 

9. A method as claimed in any one of the preceding cla ims 
wherein any unmatched elements of said first and second 
surfaces, being elements that could not be matched, are 
assigned thicknesses being the average of the thicknesses 

15 of adjacent matched elements where such adjacent matched 
elements exist f or of adjacent unmatched elements where 
such adjacent matched elements do not exist and said 
adjacent unmatched elements have been assigned thicknesses. 

20 10. A method as claimed in claim 9 wherein each element of 
an edge surface, being a Burface between said first and 
second surfaces, and adjacent to either of said first or 
second surface is assigned a thickness proportional to the 
thickness of the element of said first or second surface to 

25 which said element of said edge surface is adjacent. 

11. A me t ho d as claimed in claim 10 wherein each said 
element of an edge surface is assigned a thickness between 
0.5 and 1.5 times said thickness of the element of said 

30 first or second surface to which said element of said edge 
surface is adjacent. 

12. A method as claimed in claim 11 wherein each said 
element of an edge surface is assigned a thickness between 

35 0.7 and 0.9 times said thickness of the element of said 

first or second surface to which said element of said edge 
surface is adjacent. 
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13. A method as claimed in claim 12 wherein each said 
element of an edge surface is assigned a thickness 
approximately 0.75 times said thickness of the element of 

5 said first or second surface to which said element of said 
edge surface is adjacent. 

14. A method as claimed in claim 10 wherein each element 
of an edge surface not adjacent to said first or Second 

10 surface is assigned a thickness being the average of the 
thicknesses of adjacent elements of said edge surface that 
have been assigned thicknesses. 

15. A method as claimed in any one of the preceding claims 
15 wherein flow is simulated at a rate directly proportional 

to a desired flow rate for the object. 

16. A method as claimed in claim 15 wherein said rate is 
proportional to the ratio of computational domain volume of 

20 said object to real volume of said object. 

17. A method as claimed in claim 16 wherein said rate is 
substantially equal to the ratio of said computational 
domain volume to said real volume. 

25 

18. A method as claimed in any one of the preceding cla ims 
wherein said method is performed with first and second 
representations of said first and second surfaces 
respectively comprising first and second meshes or lattices 

30 respectively, wherein said elements are interstices of said 
first and second meshes or lattices. 

19. A method as claimed in any one of the preceding claims 
wherein said elements are triangular. 

35 

20. A method as Claimed in any one of claims 1 to 18 
wherein said elements are quadrilateral. 
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21. A method as claimed in claim 19 wherein said elements 
are substantially equilateral. 

5 22. A method as claimed in claim 18 wherein said method 
includes creating said first and second representations. 

23. A method as claimed in either claim 18 or 22 wherein 
said method includes creating improved representations of 

10 said first and second surfaces , whereby said elements are 
elements of said improved representations and said me thod 
is performed with said improved representations. 

24. A method as claimed in claim 18 wherein said first and 
15 second representations are, or are a part of, a 

representation or representations for stereolithography of 
said object. 

25. A method as claimed in any one of the preceding claims 
20 wherein said method is performed by a con^mter r unn ing a 

computer program encoding said method for simulating fluid 
flow. 

26. A method as claimed in any one of the preceding claims 
25 wherein said method includes corrections for non- isothermal 

temperature fields and/or non-Newtonian fluids. 

27. A method for simulating fluid flow within a three 
dimensional object having first and second generally 

30 opposed surfaces including: 

providing or creating first and second 
representations of said first and second surfaces 
respectively, 

creating first and second improved 
35 representations from said first and second representations 
respectively, 

matching each element of said first improved 
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representation of said first surface with an element of 
said second improved representation of said second surface 
between which a reasonable thickness may be defined, 
wherein matched elements of said first improved 
5 representation constitute a first set of matched elements 
and matched elements of said second improved representation 
constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using each set of said 
10 matched elements, whereby said injection point is linked tQ 
all locations on said first and second improved 
representations from which flow may emanate such that 
resulting flow fronts along said first and second improved 
representations are synchronized. 

15 

28. A method as claimed in claim 27 wherein said first and 
second representations are, or are a part of, a 
representation or representations for stereolithography of 
said object. 

20 

29. A method as claimed in either claim 27 or 28 wherein 
said first and second improved representations comprise 
small equilateral triangular elements. 

25 30. A method for simulating fluid flow within a three 
dimensional object having first and second generally 
opposed surfaces including: 

matching each element of said first surface with 
an element of said second surface between which a 

30 reasonable thickness may be defined, wherein matched 

elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements and 
said elements are substantially equilateral triangles, 

35 specifying a fluid injection point, 

performing a flow analysis using each set of said 
matched elements, whereby said injection point is linked to 
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all locations on said first and second surfaces from which, 
flow may emanate such that resulting flow fronts along said 
first and second surfaces are synchronized, wherein said 
first and second representations are, or are a part of, a 
5 representation or representations for stereolithography of 
said object. 

31. A method as claimed in claim 30 wherein said injection 
point is one of a plurality of injection point*. 

10 

32. A computing device provided with or running a computer 
program encoding a method for simulating fluid flow as 
claimed in any one of the preceding claims. 

15 33. a computer storage medium provided with a computer 
program embodying a method for simulating fluid flow as 
claimed in any one of claims 1 to 31. 

34. A method for modelling a three dimensional object 
20 including: 

specifying first and second generally opposed 
surfaces of said object, 

forming first and second representations Of said 
first and second surfaces respectively, wherein said first 
25 and second representations each comprise a plurality of 
elements, 

matching pairs of elements of said first and 
second surfaces between which a reasonable thickness may be 
defined. 

30 

35. A method as claimed in claim 34 wherein said first and 
secoxxd representations comprise first and second meshes or 
lattices respectively, wherein said elements are 
interstices of said first and second meshes* or lattices. 

35 

36. A method as claimed in either claim 34 or 35 wherein 
said elements are triangular. 
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37. A method as claimed in claim 36 wherein eaid elements 
cure substantially equilateral. 

5 38. A method as claimed in claim 34 wherein said elements 
are quadrilateral. 

39. A method as claimed in any one of claims 34 to 38 
wherein each element of each of said matched pairs of 

10 elements is assigned respectively said thickness. 

40. A method as claimed in claim 39 wherein unmatched 
elements of said first and second surfaces are assigned 
thicknesses being the average of the thicknesses of 

15 surrounding, matched elements of said first and second 
surfaces . 

41. a method as claimed in claim 40 wherein any unm atched 
elements of said first and second surfaces, being elements 

20 that could not be matched, are assigned thicknesses being 
the average of the thicknesses of adjacent snatched elements 
where such adjacent matched elements exist, or of adjacent 
unmatched elements where such adjacent matched el emen ts do 
not exist and said adjacent unmatched elements have be en 

25 assigned thicknesses. 

42 . A method as claimed in claim 41 wherein each element 
of an edge surface, being a surface between said first and 
second surfaces, and adjacent to either of said first or 
30 second surface is assigned a thickness proportional to the 
thickness of the element of said first or seco nd surface to 
which said element of said edge surface is adjacent. 

43. A method as claimed in claim 42 wherein each said 
35 element of an edge surface is assigned a thickness between 
0.5 and 1.5 times said thickness of the element of said 
first or second surface to which said element of said edge 
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surface is adjacent. 

44. A method as claimed in claim 43 wherein each said 
elem ent of an edge surface is assigned a thickness between 

5 0.7 and 0.9 times said thickness of the element of said 
first or second surface to which said element of said edge 
surface is adjacent. 

45. A method as claimed in claim 44 wherein each said 
10 element of an edge surface is assigned a thickness 0.75 

times said thickness of the element of said first or second 
surface to which said element of said edge surface is 
adjacent. 

15 46. A method as claimed in claim 42 wherein each element 
Cf an edge surface not adjacent to said first or second 
surface is assigned a thickness being the average of the 
thicknesses of adjacent elements of said edge surface that 
have been assigned thicknesses. 

20 

47. A method for simulating fluid flow within a three 
dimensional object having first and second generally 
opposed surfaces including: 

matching each element of said first surface with 
25 an element of said second surface between which a 

reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements, 
30 specifying a fluid injection point, 

performing a flow analysis using each set of said 
matched elements, and 

synchronizing flow fronts resulting from said 
flow analysis along said first and second surfaces. 

35 

48. A method as claimed in claim 47 wherein said flow 
fronts are synchronized from substantially the Commencement 
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of Baid flow analysis, 

49 . A method as claimed in claim 47 wherein said flow 
front* are first synchronized after the commencement of 

5 said flow analysis. 

50. A method as claimed in any one of claims 47 to 49 
wherein said injection point remains linked to all said 
locations at substantially all times in said flow analysis 

10 subsequent to being first, so linked. 

51. A method as claimed in any one of claims 47 to 50 
wherein said injection point is one of a plurality of 
injection points. 

15 

52. A method as claimed in any one of claims 47 to 51 
wherein said synchronization of said flow fronts is checked 
periodically. 

20 53. A method as claimed in claim 52 wherein said checking 
is performed at each time step. 

54. A method as claimed in any one of claims 47 to 53 
wherein said flow fronts are synchronized if found not to 
25 be or no longer to be synchronized. 

55 . A method as claimed in any one of claims 47 to 54 
wherein said first and second generally opposed surfaces 
are one of a plurality of pairs of generally opposed 
30 surfaces. 

56. A method as claimed in any one of claims 47 to 55 
wherein any unmatched elements of said first and second 
surfaces, being elements that could not be matched, are 
35 assigned thicknesses being the average of the thicknesses 
of adjacent matched elements where such adjacent matched 
elements exist, or of adjacent unmatched elements where 
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such adjacent snatched elements do not exist and said 
adjacent unmatched elements have been assigned thicknesses. 

57. A method as claimed in claim 56 wherein each element 

5 of an edge surface, being a surface between said first and 
second surfaces, and adjacent to either of said first or 
second surface is assigned a thickness proportional to the 
thickness of the element of said first or second surface to 
which said element of said edge surface is adjacent. 

10 

58. A method as claimed in claim 57 wherein each said 
element of an edge surface is assigned a thickness between 
0 . 5 and 1 . 5 times said thickness of the element of said 
first or second surface to which said element of said edge 

15 surface is adjacent. 

59. A method as claimed in claim 58 wherein each said 
element of an edge surface is assigned a thickness between 
0.7 and 0.9 times said thickness of the element of said 

20 first or second surface to which said element of said edge 
surface is adjacent. 

60. A method as claimed in claim 59 wherein preferably 
each said element of an edge surface is assigned a 

25 thickness approximately 0.75 times said thickness of the 
element of said first or second surface to which said 
element of said edge surface is adjacent. 

61. A method as claimed in claim 60 wherein each element 
30 of an edge surface not adjacent to said first or second 

surface is assigned a thickness being the average of the 
thicknesses of adjacent elements Of said edge surface that 
have been assigned thicknesses. 

35 62. A method as claimed in any one of claims 47 to 61 

wherein flow is simulated at a rate directly proportional 
to a desired flow rate for the object. 
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63. A method as claimed in claim 62 wherein said rate is 
proportional to the ratio of computational domain volume of 
said object to real volume of said object. 

5 

64. a method as claimed in claim 63 wherein said rate is 
substantially equal to the ratio of said computational 
domain volume to said real volume. 

10 65. A method as claimed in any one of claims 47 to 64 
wherein said method is performed with first and second 
representations of said first and second surfaces 
respectively comprising first and second meshes or lattices 
respectively, wherein said elements are interstices of said 

15 first and second mesheB or lattices. 

66. A method as claimed in any one of claims 47 to 65 
wherein said elements are triangular. 

20 67. A method as claimed in any one of claims 47 to 65 
wherein said elements are Quadrilateral. 

68. A method as claimed in claim 66 wherein said elements 
are substantially equilateral . 

25 

69. A method as claimed in claim 65 wherein said method 
includes creating said first and second representations. 

70. A method as claimed in claim 65 wherein said method 
30 includes creating improved representations of said first 

and second surfaces, whereby said elements are elements of 
said improved representations and said method is performed 
with said improved representations. 

35 71. A method as claimed in claim 65 wherein said first and 
second representations are, or are a part of, a 
representation or representations for stereolithography of 
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said. object* 

72. A method as claimed in any one of claims 47 to 71 
wherein said method includes corrections for non- isothermal 

5 temperature fields and/or non-Newtonian fluids. 

73. A method as claimed in any one of claims 47 to 72 
wherein said method is performed by a computer running a 
computer program encoding said method for simulating fluid 

10 flow. 

74. A method for simulating fluid flow within a three 
dimensional object having first and second generally 
opposed surfaces including: 

15 matching each element of said first surf ace with 

an element of said second surface between which a 
reasonable thickness may be defined/ wherein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
20 surface constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using said first set 
of matched elements, 

adapting said flow analysis to said second set of 
25 matched elements, and 

synchronizing flow fronts resulting from said 
flow analysis and said adaptation of said flow analysis 
along said first and second surfaces. 

30 75. A method as claimed in claim 74 wherein said method is 
performed with first and second representations of said 
first and second surfaces respectively comprising first and 
second meshes or lattices respectively, wherein *jaid 
elements are interstices of said first and second meshes or 

35 lattices. 



76. A method as claimed in either claim 74 or 75 wherein 
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said .elements are triangular. 

77. A method as claimed in either claim 74 or 75 wherein 
said elements are quadrilateral. 

5 

78. A method as claimed in claim 76 wherein said elements 
are substantially equilateral. 

79. A method as claimed in any one of claims 74 to 78 

10 wherein said method includes creating said first and second 
representations . 

80. A method as claimed in any one of claims 74 to 79 
wherein said method includes creating improved 

15 representations of said first and second surfaces, whereby 
said elements are elements of said improved representations 
and said method is performed with said improved 
representations . 

20 81. A method as claimed in any one of claims 1, 27 or 47 
wherein said synchronization comprises matching pressure 
and temperature. 



iuh 
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Figure 2 
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Figure 3 
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Figure 5 
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Figure 6 
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Figure 7a 




Figure 7b 
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Figure 9 
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Figure 12 
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Figure 13 
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Figure 17 
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